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Abstract: We study the connection between quark and lepton mass matrices in a 
supersymmetric SO(IO) GUT model in six dimensions, compactified on an orbifold. 
The physical quarks and leptons are mixtures of brane and bulk states. This leads 
to a characteristic pattern of mass matrices and high-energy CP violating phases. 
The hierarchy of up and down quark masses determines the CKM matrix and most 
charged lepton and neutrino masses and mixings. The small hierarchy of neutrino 
masses is a consequence of the mismatch of the up and down quark mass hierar- 
chies. The effective CP violating phases in the quark sector, neutrino oscillations 
and leptogenesis are unrelated. In the neutrino sector we can accomodate naturally 
sin 6123 ~ 1, sin 6*13 < 0.1 and mi < ma ~ ^/Am^ < mg ~ A/Am|^. 
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1. Introduction 

Grand unified theories (GUTs) appear to be the most promising framework [1,2] to 
address the still challenging question of quark and lepton masses and mixings. During 
the past years new results from neutrino physics have shed new light on this problem, 
and the large differences between the mass hierarchies and mixing angles of quarks, 
charged leptons and neutrinos impose strong constraints on unified extensions of 
the Standard Model (SM) [3,4]. Massive neutrinos are most easily incorporated in 
theories with right-handed neutrinos, which leads to SO(IO) as preferred GUT gauge 
group [5, 6]. 

Higher-dimensional theories offer new possibilities to describe gauge symmetry 
breaking, the notorious doublet-triplet splitting and also fermion masses. A simple 
and elegant scheme is provided by orbifold compactifications which have recently 
been considered for GUT models in five and six dimensions [7-12]. In this paper 
we analyse in detail the connection between quark and lepton mass matrices in the 
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six-dimensional (6D) GUT model suggested in [13], for which also proton decay [14], 
supersymmetry breaking [15] and gauge coupling unification [16] have been studied. 
An alternative SO(IO) model in five and six dimensions has previously been studied 
in [17]. For a recent discussion of CP violation in a 5D orbifold GUT model, see [18]. 

An important ingredient of orbifold GUTs is the presence of split bulk multiplets 
whose mixings with complete GUT multiplets, localised at the fixed points, can sig- 
nificantly modify ordinary GUT mass relations. This extends the known mechanism 
of mixing with vectorlike multiplets [19-21]. Such models have a large mixing of 
left-handed leptons and right-handed down quarks, while small mixings of the left- 
handed down quarks. In this way large mixings in the leptonic charged current are 
naturally reconciled with small CKM mixings in the quark current. 

Our model of quark and lepton masses and mixings relates different orders of 
magnitude whereas factors C(l) remain undetermined. Hence, we can only discuss 
qualitative features of quark and lepton mass matrices. Recently, orbifold compact- 
ifications of the heterotic string have been constructed which can account for the 
standard model in four dimensions and which have a six-dimensional GUT structure 
as intermediate step very similar to familiar orbifold GUT models [22-24]. In such 
models the currently unknown C(l) factors are in principle calculable, which would 
then allow for quantitative predictions. 

The goal of the present paper is twofold: As a typical example, we first study 
the model [13] in more detail and explicitly compute the mass eigenstates, masses 
and mixing angles. Second, we investigate the question of CP violation, both in 
the quark and lepton sector and possible connections between the two. In previous 
studies, CP violation has mostly been neglected assuming that, barring fortunate 
cancellations, the phases and mixings are practically independent. Nevertheless this 
question and the fiavour structure are strongly interconnected, and we will see that a 
specific pattern of mass matrices can give a distinct signature also in the CP violation 
invariants. 

This paper is organised as follows: In Section [2[ we describe the 6D orbifold 
GUT model and the diagonalisation of the mass matrices defining the low energy 
SM fermions. In Section [^ we discuss the CP violation in the quark sector, whereas 
Section |[ is devoted to the CP violation in the leptonic sector. Conclusions are 
given in Section [^. Two appendices provide details to the computation of the mass 
eigenstates and CP violation in extensions of the SM. 

2. SO(IO) Unification in six dimensions 

We study an SO(IO) GUT model in 6D with = 1 supersymmetry compactified 
on the orbifold TV(Z2 x x ^2^) [11,12]. The theory has four fixed points, Oj, 
Ops, Oqq and Og, located at the four corners of a 'pillow' corresponding to the two 
compact dimensions (cf. Fig. 0). The extended supersymmetry is broken at all fixed 
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Figure 1: The three SO(IO) subgroups at the corresponding fixed points (branes) of the 
orbifold TV(Z^ X x 'Z^^). 



points; in addition, the gauge group SO(IO) is broken to its three subgroups Gps = 
SU(4) X SU(2) X SU(2); Ggg = SU(5) x U(l)^; and flipped SU(5), Gfl = SU(5)' x 
U(l)', at Ops, Ogg and Og, respectively. The intersection of all these GUT groups 
yields the standard model group with an additional U(l) factor, GgM' = SU(3) x 
SU(2) X U(l)y X U(l)y,, as unbroken gauge symmetry below the compactiflcation 
scale. 

The field content of the theory is strongly constrained by imposing the cancella- 
tion of irreducible bulk and brane anomalies [25]. The model proposed in Ref. [13] 
contains three spinors ipi{16), i = 1 ... 3, as brane fields as well as six vectorial fields 
Hj{10), j = 1...Q, and two pairs of spinors, $(16) + $''(16) and <^(16) + 0''(16) as 
bulk hypermultiplets. 

The massless zero modes A^($) and A^''(<l>'') acquire vacuum expectation values 
(vevs), vn = {N) = {N'^), breaking B — L and thus Gsm' to Gsm- The breaking scale 
is close to the compactification scale so that v%/M^, ~ 10^^ GeV, where M^. is the 
cutoff of the 6D theory. At the weak scale, the doublets Hd{Hi) and Hu{H2) acquire 
vevs, Vi = {Hd) and V2 = {Hy)^ breaking the electroweak symmetry. 

The three sequential 16-plets are located on the three branes where SO(IO) is 
broken to its three GUT subgroups; in particular, we place ipi at Ogg, "^2 at Og and 
'03 at Ops- The parities of H^, Hq, cj), and (p^ are chosen such that their zero modes, 

m = (l^^) , = (l^f) , dim , d,m , (2.1) 

have the quantum numbers of a lepton doublet and antidoublet as well as anti-down 
and down-quark singlets, respectively. Both L(0) and L^^cj)'^) are SU(2)^ doublets. 
Together these zero modes act as a fourth vectorial generation of down quarks and 
leptons. 

The three 'families' ipi are separated by distances large compared to the cutoff 
scale M^,. Hence, they can only have diagonal Yukawa couplings with the bulk Higgs 
fields; direct mixings are exponentially suppressed. The brane fields, however, can 
mix with the bulk zero modes for which we expect no suppression. These mixings 
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take place only among left-handed leptons and right-handed down quarks, leading 
to a characteristic pattern of mass matrices [13, 14]. 
The mass terms assume the characteristic form, 

W = Uimyi + d^m^d'^p + elmlpep + vlm^pVp + \ vtmfvt , (2-2) 

where latin indices only span 1, 2, 3, while greak indices include the forth generation 
states. The up quark and Majorana neutrino mass matrices, m" and , are diagonal 
3x3 matrices, 



'h'l-^V2 

m'^ = I /i22f2 

/i^gVsy 



m 



N 



fuN 4l 
I ""11 M, 



V 





'^22 M* 








hg^^J 



(2.3a) 



Since is part of an SU(2)^ doublet, it cannot couple to the other SM singlets in 
■01 via the B — L breaking field. Furthermore, there is no other coupling giving it a 
direct Majorana mass. 

The Dirac mass matrices of down quarks, charged leptons and neutrinos, m*^, rrf 
and m^, respectively, are 4x4 matrices instead, due to the mixing with the bulk 
field zero modes, 



m 



m 



(Kivi gij^v,\ 
hi^v, gij^v^ 

\flVN f2VN hVN M'^ ) 



fhf^vi hl4Vi\ 
hl2Vi h24Vi 
h^^vi hl^vi 
\ M[ Mi, Ml M[ j 



(2.3b) 



(h^^V2 h^^V2\ 

/l22'^2 h2^V2 
hl^V2 hgv2 

\ Mi Mi, Ml Mi J 



(2.3c) 



up to corrections 0{vlf/M^). The diagonal elements satisfy four GUT relations 
which correspond only to the local unbroken groups, i.e., SU(o), fiipped SU(o) and 
Pati-Salam subgroups of SO(IO). The hypothesis of a universal strength of Yukawa 
couplings at each fixpoint leads to the identification of the diagonal and off-diagonal 
elements of m"/ tan/?, m'^, rrf, and m^/ tan/3, where tan/? = ^2/^1, up to coeffi- 
cients of order one. This implies an approximate top-bottom unification with large 
tan 13 and a parametrisation of quark and lepton mass hierarchies in terms of the six 
parameters /ij and /Ij. 

The crucial feature of the matrices mf', rrf and are the mixings between the 
six brane states and the two bulk states. The first three rows of the matrices are 
proportional to the electroweak scale. The corresponding Yukawa couplings have to 
be hierarchical in order to obtain a realistic spectrum of quark and lepton masses. 
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This corresponds to different strengtlis of the Yukawa couplings at the different fixed 
points of the orbifold. The fourth row, proportional to M'^, M' and Vn, is of order 
the unification scale and, we assume, non-hierarchical. 
The mass matrices m'^. 



and are of the common form 
/ /ii /ii \ 

/i2 /i2 

J)^ ^ ^ 
\Mi M2 M3 M^J 



m 



(2.4) 



where = C(fi,2) 

the unitary matrices 



and Mi = (9(Mgut)- This matrix can be diagonalised using 



m 



(2.5) 



where the matrices t/4 and V4 single out the heavy mass eigenstate, that can then 
be integrated away, while U3, and V3 act only on the SM flavour indices and perform 
the final diagonalisation also in the 3x3 subspace. The explicit expressions for the 
mixing matrices and the mass eigenstates are given in Appendix 0. 

The parameters in the matrix Eq. ( |2.4| ) are generally complex; however, we can 
absorb seven phases with appropriate field redefinitions and choose the remaining 
three physical phases to be contained into the diagonal parameters /Xj, 



m 



1^: 



2 e 



iB2 



\ 






Ms 



7^2 



(2.6) 



M4/ 







Ml M2 

This is the maximal number of physical phases for four generations of Dirac fermions, 
given as usual by — l)(n — 2)/2 for n generations, so our texture above does not 
reduce the CP violation from the typical n = 4 case. We will see that the phases 
survive in the low energy parameters, but that only one combination defines the 
single phase characteristic of three generations. 

With this choice, the matrix V4 is real, while t/4 contains complex parameters; 
however, the imaginary part is suppressed by /M so that their effect on the low 
energy CP violation is negligible as long as the mass of the heavy eigenstate is large 
compared to the electroweak scale. From the unification of the gauge couplings, we 
expect indeed M to be of the order of the GUT scale [16]. Then the discussion of 
the low energy CP violation, which would in general be characterised by many CP 
invariants [26,27], reduces to the case of three light generations (see Appendix P). 

The effective mass matrix is given by m, the 3x3 part of 

^/il(V4)lj +Jll{Vi)4i 



m 



UlmVi 




VM2 



m = I fi2iV4)2j +/i2(^4)4i 



(2.7) 
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m 



in terms of the parameters in Eq. (|2^) , it reads 

Ui — / = — til — , = u — P'l^^ — / — Ui^^ — / - 

^^yM^fMf ''^'^ ^/Iq+ul My/Ml+Ml Ma/m2+M| 

\ Y^M2+M| y/Mf+Mi My/Mi+Mi ^^Ma/m^TmI / 

As any matrix, m can be transformed into upper triangular form just by basis redef- 
inition on the right, 

rn = m Vs = Ji2 0/^2 • (2.8) 
V /is / 

This form is particularly suitable in the case of the down quarks, where V3 acts on 
the right-handed quarks and disappears from the low energy Lagrangian due to the 
absence of right-handed current interactions. Note that we can reshuffle the phases, 
reabsorbing three of them into the unitary transformation V3, but we are still left 
with three complex parameters. We can exploit this freedom to obtain real diagonal 
elements /22, /Is and 7/2i, while a, (3, and /2i remain complex. 

On the other hand, we can still redefine two phases on the left-hand side, keeping 
an overall phase free, with a diagonal matrix 



L3 



diag (e-^^^ e-*^^ 1) . (2.9) 



This transformation allows us to shift the phase of JLi into 7, which will be convenient 
later in the limit where 7 vanishes. Again, such a phase shift does not reduce the 
number of complex parameters in the down quark matrix, which remains three. 
Moreover, this reparametrisation does not change the CKM matrix, since the up 
quark mass matrix is diagonal and so such phase transformation can be compensated 
by an identical one for both Ui and u^. 

The matrix V3 differs from the upper 3x3 part of the diagonalising matrix 
V3 = VJjVJj'; however, they are very similar in the hierarchical case. The relation 
between these two can be found in Appendix together with the general expression 
for [/s, the 3x3 part of which is the CKM matrix. 

For the leptons, it is the matrix V4V3 that acts on the left-handed states, so the 
mismatch between the charged leptons and neutrinos (see Eq. (p.3c|) ) basis appears in 
the charged current interaction and the definition of the flavour neutrino eigenstates. 
However, the matrix V4 which contains large mixing angles and rotates away the 
heavy eigenstate is the same for charged leptons and neutrinos since the heavy state 
is an SU(2)^ doublet. Therefore the PMNS matrix will be given only by the mismatch 
between the V3 ~ V3 matrices for charged leptons and neutrinos. 
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The complete expressions for the parameters in fn are given in Appendix ^ in 
this section, we will only consider the limit of small /xi as well as small Jli and/or 
/i2- For Hi = Jli = 0, the first row simply vanishes, whereas for fii = = 0, the 
two first rows of the mass matrix are aligned (see Eq. (|2?6|)). Therefore both cases 
correspond to vanishing down-quark and electron mass. 

Since JI1/JI2 gives Vus, we focus on the case /ii = /i2 = 0, where^ 





= (3 = 


/f2 / 


/is 


a 














El 


- El 


_ 1/^3 




Jli 




JJ'3 





M4 /i^Mg + /I3M4 



M 



Ml + M| 



M 



The eigenvalues of the heavier states are given by 



7 = 



(2.10) 



m 



2 7t2 



1^2 



where ^ ~ V^^ 

/i2 



(2.11a) 
(2.11b) 



In this limit, only one single physical CP violating phase survives, even in the 4x4 
picture; it is contained in /is and so in a and (3 (see Eq. (|2.10|) ). We will see, however, 
that this single phase is not sufficient to have low-energy CP violation. 

The down-quark mass is indeed very small, so we will use these expressions as 
the order zero approximation, together with the corrections proportional to |/i2| //22, 
which determine the masses of the down-quark and the electron. Our expansion 
parameter will therefore be of the order of the mass ratio of the down and strange- 
quark, m^/ms. In fact, for <^ \jj,2\ we have at leading order 



so our expansion parameter is 



fj'i— — — 

fJ'2 /Lis 



/i2 



Vus |/Lt2| 



(2.11c) 



0.23 



(2.12) 



The mass ratio of electron and muon is much smaller than the ratio of down 
and strange quark. This implies (/i2/ii//i2)e (/i2/Ui//I|)d. Assuming that the 
difference is due to the smallest matrix elements, this indicates (/i2)e/(/i2)d ^ 1 
and/or (/^i)e/(/^i)d ^ 1 for (/i2)e — ijl2)d- This fact can easily be accommodated, 
as we see in Eqs. ( p.3|) : the presence of the second generation on the fiipped SU(5) 



^As mentioned above, it is instructive to choose the basis in which Jii is real and the vanishing 
parameter 7^1 complex. Then it is obvious that we are left with only two complex parameters in 
TO, namely a and /3, containing the same phase. 
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brane leads to different values of /i2 for the down quarks and charged leptons and 
the parameter Jli stems from different couplings in the superpotential. 

While we derived the fermion mass matrices ( p.3| ) within a specific model, they 
can also arise in other models, where additional matter is present at the GUT (or 
compactification) scale. Thus we could take these matrices as a starting point for 
the following discussion, leaving open the question of their origin. 



3. CP violation in the quark sector 

We will first consider the CP violation in the quark sector. As we have seen in the 
previous section, our effective 3x3 down quark mass matrix contains three phases 
as a remnant of the original 4x4 matrix, with the dominant complex element being 
a/22. We will now derive the combination of the three phases, which plays the role 
of the CKM phase. 

To describe CP violation for three generations, as is the case in the SM, it is 
convenient to use the Jarlskog invariant [28], Jg, which is given by 

6tA^^A^^Jg = tT[H^,H,f = Q \mti [HlHjH^H,), (3.1) 

where H = mm) and 

= [ml - ml) (ml - mj) (ml - mj) ; (3.2) 

note that has mass-dimension six. In our model, the up quark mass matrix is 

diagonal, as is Hu. Then the invariant strongly simplifies and reads 

rl2 0-23 o'31^ 



ImiH^'HlHf) 



' (3.3) 



It is clear from this expression, that any diagonal phase transformation of m on the 
left does not have any effect on the Jarlskog invariant. 

As discussed in Appendix |B|, we can use the effective 3x3 mass matrix 
^eff _ = f7|f7|t By means of Eq. (|2.8|) , we obtain 

(\jLxf (l + + Yif) JI1JI2 (1 + iui/i3/3\ 
/I1/22 (1 + a/5*) Jll{l + \af) /22/23a , (3.4) 
Kf^sf^* /22/i3a* Jll J 

where JI2 and JI3 are real parameters, as displayed in Eqs. (|A.4|) . Then we have 



Im 



(Hf)'" (Hf)'' (Hf)'' = \jl,\'jlljlllmaP* (1 + a*/5) (3.5) 



Itt |2 -2^2 ■ 



\fiif fJ'lfil Ima/3* 
JI2J1I Im[(a/22) {Pjli)*Jli] 
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We see that the Jarlskog invariant is always independent of the argument of 7 and it 
vanishes in the hmit /ii, /Ii — > such that jui = 0. As we might expect, Jq vanishes 
for a = /9 as well, i.e., in the limit /ii, ^2 ~^ 0. 

So the presence of a single phase in a is not sufficient to give CP violation in the 
low energy: this phase cancels out in the Jarlskog invariant. This effect stems from 
the alignment of the vectors in flavour space; however, even in the case of vanishing 
flrst generation mass, the corresponding eigenvector does not decouple from the other 
two and the mixing matrix does not reduce to the two-generational case. In fact, the 
CKM matrix is given by (see Appendix 



CKM 



(m, = 0) 



v 



1 




/ilQ\ 




J12 


/is 




1 


/i2Q 


/i2 




/i3 





/i3 


1 ) 



n 



(3.6) 



Hence, we cannot conclude that the CP effects disappear due to the reduction of 
the system to two generations, nor to the mass degeneracy between quarks. Instead 
the absence of low energy CP violation is caused by the particular texture of m in 
exactly the same basis for the left-handed quark doublet, where the up quark matrix 
is diagonal. This feature is similar to the absence of CP violation in 4D SO(IO) 
constructions, where a single ten-dimensional Higgs field generates fermionic masses, 
yielding a trivial CKM matrix. Note that there is still some CP violation effect 
arising from the dominant phase ^3 in /^s, but it is only apparent in the mixings 
involving the fourth heavy state. 

Now, the down quark is not massless and the real physical case corresponds to 
non-zero /ii, /i2 and yUi. From the up quark phenomenology, we know that /ii : yU2 
is similar to the mass ratio of up and charm-quark [13]; in addition, : /X2 is fixed 
by the Cabibbo angle. We will therefore focus on the linear terms in /i2 and keep 
/ii ~ 0. 

As is apparent in Eq. (|3.5|) , contributions to Jq come from the complex quantities 
a/i2, /9/ii, and jui; however, is independent of /i2 (see Eq. (|A.4]) ), 



= 

The first order terms are 

5(a/22' 



/is 
/is 



M4 /I3M4 + /i*M3 



M 



(3.7) 



M2 /i3M4 + /i*M3 
-/i2^ ^ , 

M JI3M 



= Hi 



/i* /i3 M2 fisMs - /i^M4 



/i2/i3 M 



^We can exploit the phase transformation P^l ( |2.9| ) to absorb the phases of ^1, a and make all 
elements of the CKM matrix real showing explicitly that the CP violation disappears. 
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and the Jarlskog invariant reads 



Ja 



1 - 



Mi 
M2 



im ^ ^ + 



M2 



Im 



/i2 



(3.9) 



We see that Jg vanishes if either ^2 or /is vanish, so two complex quantities are 
needed to obtain CP violation at low energies. 

It is instructive to calculate H^^ also from the matrix m, Eq. (|2.7] ). Here we notice 
that the off-diagonal elements of such matrix are relatively simple since we can exploit 
the unitarity of the matrix V4, which gives YM:=i(^4)ik(^4:)% = ^ij ~ (^4)i4(^)j4- So 
we have for i j 



(Hf) = fiifij (1 - aia*) 



(3.10) 



from which we get the simple expression 

, 12 / TTom 23 / TTPfTX 31 



Im 



{Hfr {Hfr {HfY 



~2~2~2 



E (1 

cycl. perm ijk 



) Im {a*ak) . (3.11) 
0, the 



In the limit of vanishing /ij, we see that Oj = M4/M; thus for /ii = /i2 
expression simplifies to 



Im 



{HffiHfr (Hf) 



23 



,31 



~2~2~2 I 1 I 
/^l/^2^3 1 + 



M2 



For /ii = but /i2 7^ 0, we then obtain 

12 / TTofF\ 23 / TTpfFN 31 



Im 



{Hfr {Hfr {HfY 



P1P2A-3 ^2 



Im 



M2 



,M4 , M4 

M M 







Im 



M3M4 |Ai3| 



Im 



/i2 
/i2 



V/^3/^2/ 

M2M4 |/i2| 



(3.12) 



M2 ^^2 



Im 



/is 
/W3 



M2 /is 

The complete expression for Jq is displayed in Eq. ( |A.19|) ; the dominant terms are 
exactly those given in Eq. (|3.9|). 

For degenerate heavy masses M, the result simplifies to 

1 /X?/i2 



J, 



" 16 A^2 



3 /i3 Im (/ig/i^) + I /i3 1 Im(/i2) 



(3.13) 



Note that the numerical factor, ^, is minimal for degenerate M. Due to the hierarchy 
of the down quarks, ~ mimf ~ Jill^t- So we finally obtain, substituting the 

order of magnitude of the parameters, with /I3 ~ IfJ'sl, 



mdrris i 
/ ~ 1/ 

'Jq — ^ us 9 , , — 

ml 4v/2 



(3 sin (03 - O2) + sin02) ^ 10"' (3 sin (^3 - ^2) + sin02) • 



(3.14) 
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This is the right order of magnitude; the current experimental value is = 3 x 10^ 
[29]. From Eq. (|3.14|) we can conclude that a single complex parameter, with the 
other two vanishing, is not enough to have low-energy CP violation in the quark sector 
and that the CKM phase is a combination of the high-energy phases 6i weighted by 
mass hierarchies. Moreover, maximal phases seem to be needed to give the large 
low-energy phase observed. 



4. CP violation in the leptonic sector 

The charged lepton and Dirac neutrino mass matrices can be transformed like the 
down quark mass matrix. The heavy state is an SU(2)j^ doublet, so V4 singles out 
the same state for charged leptons and neutrinos. 
The effective 3 x 3-matrices read (cf. Eq. (|277|)) 

fl2{Vi)2] + MyA)ij y = P2{Vi)2j + p2{Vi)ij ■ (4.1) 

Pziyi)?,3 + V'?.iyi)Aj/ \P3(V'4)3i + P3(V'4)4i/ 

Within our model we assume the hierarchical patterns of pi and as well as pi and 
Pi {i = 1..3) to be the same as for down quarks. The precise values, however, can 
be different since they originate from different Yukawa couplings, see Eqs. (|2.3c|) . 
Again, we choose the couplings between the brane states, /i, and Pi, complex. 

Although some of the charged lepton and down quark parameters, namely Hi 
and ps, are related by GUT symmetries, the corresponding phases after the redef- 
inition leading to Eq. ( |2.6| ) are completely uncorrelated. Thus, there is no direct 
relation between the CP violation in the leptonic and in the hadronic observables, 
even though, barring cancellations, we expect the leptonic CP violation to be large 
as well. Furthermore, we will see that different combinations of the phases determine 
the experimental observables. Thus even if there were relations between the phases 
in the quark and lepton sector, these would not be observable. Some correlations, 
however, could survive between charged and neutral leptons. As in the quark sector, 
we expect similar suppression for the CP violation due to the specific mass texture 
in our model. 

The discussion of the charged lepton masses closely follows the discussion of the 
down quarks in the previous section. The parameters are chosen such that they 
match the observed hierarchy, as described in Appendix |A.1|. The light neutrino 



masses, however, result from the seesaw mechanism, since we have heavy Majorana 
masses for the right-handed neutrinos. This Majorana matrix is diagonal, but can 
have complex entries (cf. Eq. (|2.3a|) ), 

'Mie^^-^i \ /Mie^^^'^is 

= I M2e^"^' = e^''^^ Mse^^^^^^s g | , (4.2) 

M^e^'^^j \ M3 
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where /S.<pij = (pi — (pj. Altogether, we have nine independent phases in the lep- 
ton sector; in the limit of small /ii and pi, they reduce to seven. Since neutrinos 
are Majorana, we have less freedom in the phase reshuffling. However, except for 
electroweak breaking effects in [4, the heavy state is effectively an SU(2)-doublet of 
Dirac fermions. This allows us to absorb some phases in the Dirac mass matrix and 
reduce the system to three generations for both charged and neutral leptons at the 
same time. In the following, we will neglect any effect of this heavy fourth generation 
doublet and concentrate on the three light generations including the right-handed 
neutrinos. We expect this approximation to be valid as long as M ~ Mqut is much 
larger than the Majorana masses Mj [16]. 

4.1 Seesaw Mechanism and EfFective Mass Matrix 

In the case of the leptons, neither m'^ nor fh^ is diagonal and therefore we will change 
the basis in order to simplify the discussion of the CP violation. Luckily, the large 
rotations of type V3, which bring the Dirac matrices into triangular form, are similar 
for charged leptons and neutrinos, thanks to the same hierarchical structure. 

To distinguish the flavour of the light neutrinos, we first act on the neutrino 
Dirac mass matrix with exactly the same V3 that transforms the charged lepton 
mass matrix into the upper triangular form, see Eq. (p.8|), and obtain 



/ Api Dpi pi\ 
= Bp2 Ep2 P2 ■ (4.3) 
\Cp3 Fps psj 

At this stage the charged lepton mass matrix is not yet diagonal, but not very far 
from it: the complete diagonalisation can be obtained by applying another nearly 
diagonal rotation matrix on the right, corresponding to the mismatch between V3 
and V3, and a CKM-like rotation U3, on the left as described in Appendix 0. Note 
that such a transformation from the left, as U^, in this case acts on the right-handed 
fields and leaves both H = fn'^ fn and the light neutrino Majorana mass matrix, 



eff 



(m^)^(m^)"'m^, (4.4) 



unchanged. In fact f/4 acts in very good approximation as the unity matrix on 
up to terms C(f^/M^), while f/3 just cancels out. 

So apart for the small rotation on the right needed to diagonalise H, which 
affects the CP violation in the neutrino oscillation only weakly (see Section |47^), the 



neutrino masses and mixings can be obtained from Eq. ([4.4|) , in the form 



C^g3 + B'^Q2 + A^Qi CFgs + BEq2 + ADqi C qz + Bq2 + Aq 



1 

2^ I Z?2^ I n2. 



mis = - [ CFqs + BEq2 + ADqi F'g^ + E'q2 + D'qi Fg^ + Eq2 + Dq^ 
Cqs + Bg2 + Aqi Fqs + Eq2 + Dqi Qz + Q2 + Qi 



(4.5) 
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where Qi = e '''^pf / M^. Note that the determinant of the (23)-submatrix of m^g 
is not of order instead it reads ^3^2 {F — E)^ + q^qi [F — D)^ + Q2Q1 {E — D)^, 
allowing a large solar mixing angle [30]. 

The leading part of the light neutrino mass matrix is obtained in the limit 
fii, pi ^ 0. From the general expressions (|A.24D one obtains 



1 1 



Pi = Pi- 

P2 



1 1 

' P2 



P3 M2 
1 1 



P3 = - 



P3 M2 



Ji^Mt^^ - plM^Mi -P2M2 p^M^ + pIM- 



Pa 



-P3 



p^M^M^ - pIM\ 



124 



(4.6) 



M2 + M2 + M2. 



where we have introduced M^p-y = 

In our model, the Dirac neutrino mass matrix has a hierarchical structure similar 
to the one of down quarks and charged leptons. The three smallest elements, however, 
have a considerable uncertainty. Since rUe 7^ ma, these elements cannot be equal for 
the three matrices. Inspection of m'^ suggests for pi the range between and mgVus', 
the difference is a factor 0{1). In the following we shall consider the case of small 
pi. For large pi so that \gi\ > {g^l, in the following discussion we should interchange 
Ps, Qs with pi, Qi and consider it as the dominant scale. 

We here assume ps : p2 : pi ~ ps : p2 : pi ~ m;, : : m^, which yields [13] 



pfMs 
^ piM2 



rrit 

~ — ~ 0.2 



\Qi\ 



P1M3 
Pi Ml 



ml rriu 



0.2 



(4.7) 



such that f?i ~ ^2 < f?3- Hence, in this model, the weak hierarchy in the neutrino 
sector can be traced back to the nearly perfect compensation between down and up 
quark hierarchies. 

The relation ^1 ~ Q2 implies for the two small neutrino masses \mi\ ~ |m2| 
barring cancellations or small parameters. As computed in Appendix |^, the masses 
at leading order assuming to dominate are given by 



\m2'mi\ 



-0, {l + \Ff + \Cf), 

\{F - E){A - B) + {D - E){B - C)\^ 
l + \Ff + \Cf 



(4.8) 



The light neutrino mass spectrum has normal hierarchy, and the ratio m^/ml can 
be identified with Am^^^/ Arn^^^, which is indeed consistent with observations within 
the theoretical uncertainties. 



13 



The coefficients A 
/ii, pi 0, 

A = Pi P2 /i3 Ml 
pi P2P3M 

P2P2 - P2P2 /is Ml 
P2JI2 P3 M 

P3P3 - P3P3 P2 Ml 
P2M 



F of tlie neutrino mass matrix rrtgg become in the limit 



B 



C 



D 



P3P'3 
Pi 1 1 



Pl P2 pI M2 



E = D + 



/i2 l/isi M^2 + plp-iM2 ( P3MZ 

P1P2 1 1 1 



Pl P2 P2 P3 M2 
1111 



]12P;M2M3 + fi; ( Jl3M't^ - pfMsM^ 



P3 P2 pI M2 



(psPs - PsPa) |^P2P3^i2 + P2^2 (^PsMa - pIM^ 



(4.9) 



Note that B, C, F vanish in the limiting case of equal hierarchy in the neutrino 
and charged lepton Dirac mass matrix, i.e., for Pi/pi = Pi/pi, and A is in this case 
proportional to 7/2i. In fact, if the neutrino and charged lepton vectors are perfectly 
aligned in flavour space the neutrino Dirac matrix becomes triangular at the same 
time as the charged lepton one and we cannot reproduce large neutrino mixing. 
There is though no reason to expect such alignment since the parameters pi, pi are 
not related by any GUT relation, as can be seen in Eq. ( p.3c|) . So the large neutrino 
mixing angles are not generated simply by the large LH rotation contained in the 
charged lepton's V3, but from its misalignment with the neutrinos. 

Using the relations between pi, pi and pi, and pi, pi and pi due to the hierarchical 
structure of the mass matrices in our model, one obtains the simple expressions. 



A ~ C ~ 3^ 

P2 



R P2 P2 
^ ~ — — 

P2 P2 



D ^ E 



(4.10) 



The mixing angles are computed in Appendix |A.2| ; in the case the parameters 
A, C are small, they are given by 



tan 6^23 ^ 
tan ^12 ~ 



151 



sin 



\E-F\ 
C 



13 



1 + 



+ 



\F\\ 

B{EF+l) 



\Q2\ 



2n3/2 



^'31 



(4.11) 



The atmospheric mixing angle 6*23 is naturally large; the current best fit [29, 31] 
restricts the parameter F as 0.7 < |F| < 1.4 to have it maximal. Note that F > 0.7 
can naturally be obtained even for Ipal /ps ~ {p^l /p^, as discussed in Appendix |A.2 
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For (/i2//^2)e ~ {fJ^2/ji2)d ~ 0.1 One then obtains \C\ ~ 0.1 and a value for 6^13 
close to the current upper bound. In this case though, p^l has to be suppressed with 
respect to the down quark case in order to give a consistently small me. The large 
solar mixing 6*12 can then be achieved for 5 ~ 0.1 — 1 with moderate tuning oi E — F. 

Another possibility is that a very small /i2 is called for to explain the smallness 
of the electron mass. In this case, we have naturally \A\ , \C\ ~ 0.01 and the reactor 



angle is dominated by the second term in Eq. ( |4.11|) . Then the angles 6^12 and 6*13 
depend on the same parameter B, but for the second one there is a suppression by 
02/ Qs- So in the case of hierarchical p,, both a large and small angle can be explained 
even with relatively large B. Such value for B is not unnatural, even for small /i2, 
if we accept p2 > {fJ'2)e- In this case we have sin ^13 < 0.1 correlated with the mass 
eigenvalues mi < m2 < rris. Note that in general, if all parameters A, B, and C are 
smaller than one, we obtain the prediction mi < m2, while for I? ~ 1 the two lowest 
eigenvalues are nearly degenerate. 

The largest of the heavy neutrino masses is given by M3 ~ mf / a/ Aml^^ ~ 
10^^ GeV. For the lightest heavy Majorana state the model provides the rough 
estimate Mi ~ M^m^/mt ~ 10^° GeV. 

4.2 Neutrinoless Double Beta Decay (Oz//3/3) 

The simultaneous decay of two neutrons may result in neutrinoless double beta decay, 
e.g., ''^Ge — > ''^Se + 2e. This process is currently most promising to prove the 
Majorana nature of neutrinos. The decay width can be expressed as 

r = |mee|% (4.12) 

where G is a phase space factor, the nuclear Oi//3/3 matrix element, and m^e is the 
(ll)-element of the light neutrino mass matrix. 

Since the electron mass is very small, the charged lepton mass matrix in trian- 
gular form has nearly a vanishing first row. Then the left-handed electron is already 
singled out; the remaining rotation mostly affects the (23)-block. Therefore we can 
already make an estimate of m^^ from the effective neutrino Majorana matrix, ^gg-. 
From Eq. (|4.5|) , we read off 

\mee\ = \C^Q3 + B^Q2 + Qi\ , (4.13) 

where the last term can be negleted. This result has the same form as the standard 
formula in the case of hierarchical neutrinos [32], 



m. 



v/A^sin^ ^i3e*(««-«^) + \ Ami, sin=^ cos^ ^13 



(4.14) 



where ^3 and ^2 are the two Majorana phases in the conventional parametrization of 
neutrino mass matrix ( A.34|) . 
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We can estimate the size of |mee| in our model using 

\p3\-p3, \P3\-P3, ^ ~ 1, l^sl ^ y^Aml^^, 1^2! ^ V A"^soi ' (4-15) 

P2 ^ 



which gives 



\fTlp 



P2 P2 



(4.16) 

Clearly, the last term dominates, yielding the familiar result for hierarchical neutrinos 
\mee\ < A/Am^ ~ 0.01 eV if fj,2/p2 < P2/P2. 

4.3 CP Violation in Neutrino Oscillations 

Leptonic CP violation at low energies can be detected via neutrino oscillations, which 
are sensitive to the Dirac phase of the light neutrino mass matrix. For a diagonal 
charged lepton mass matrix, the strength of Dirac-type CP violation is obtained from 
the invariant [27] 

tr [/^^ h']^ = Qt ^^}^u^i2 ^^.^23 ^^.^sij ^ ^^ -^^^ 

where h'^ = (m*^)^ m'^ and is the product of the mass squared differences of the 

charged leptons, cf. Eq. ( p.2|) . This quantity is connected to the leptonic equivalent 
of the Jarlskog invariant through 

J, = --^^lm[{hn''{hT{hT], (4.18) 

where 

^ ^^2 _ ^2^ j^^2 _ ^2^ ^^2 _ ^2^ ^ 6ml^^6ml^^ ~ |^2Hfe|^ (4.19) 

is now the product of the light neutrino mass squared differences. In the standard 
parametrisation given in Eq. ([A.34| ), 

Je = Im[(K)ii(K)22(K)l2(K);i] = I cos sin 2^13 sin 2^12 sin 2^23 sin (5 , (4.20) 

where 6 is the CP violating Dirac phase in the SM with massive neutrinos. 

The expressions (|4.18|) and (|4.20|) assume that the charged lepton mass matrix 
is diagonal. In our case, this matrix is nearly diagonal after the V3 rotation, as 
the electron mass is very small; in fact, the remaining rotation deviates from a 
unit matrix only in the 23 sector and at order O {m?^/m^^ <^ 1 (see Eq. (|A.13|) ). 
Therefore up to corrections of this order, we can use Eq. ([4.18| ) with the effective 
neutrino mass matrix mgg- given in Eq. ( |4.5| ), i.e., 

Ji = -3^ Im m^Y' {K^r iK^t] , (4.21) 
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with now /igg = (mgg)^mgg. We compute the first few terms and obtain 



.2 C*F |m3 



2 



(^eff) = , ,^,2 r'|2 + C*E {1 + F*E+C*B) Q2Ql + B*F {l + F*E+C*Br q^q, , 
{hlsf^^ \ [T! ,2 + {1+F*E+C*B) Q2qI + E* {1+F*E+C*B)* g*gs , 



F* 


|m3 


2 


1+\F 
C 1 




2 


1 + |F|^ + |C|^ 



iKsf'- , ,7J,2"''L,2 + ^ {1+F*E+C*B) Q2Ql + C (1+F*£;+C*B)* qIq^ 

(4.22) 



The leading contribution in the cyclic product of /igg-, which is oc jmsl^, is real and 
does not contribute to J^; that is to be expected since it corresponds to the limit of 
two massless neutrinos where no physical Dirac phase can be defined. In general the 
first non-trivial terms are of order |i?3|'' |^?2|^, as A^^, so that we expect \Ji\ < 1. 
We obtain in fact 

J, ~ A U ' ' > ^ I ' ' ^ Im [C*F{F - EfiB - C)] . 

(1 + |F|'+|C|') 

(4.23) 

Note that the imaginary part vanishes for E — F or B — C, when the flavour eigen- 
vectors are partially aligned. Furthermore, the contribution disappears for C = 0, so 
it is suppressed by the small reactor angle as expected. Due to the unknown param- 
eters 0{1), no useful upper bound on Ji can be derived in the general case, but we 
see that the Dirac CP phase is given by a combination of the phases of the neutrino 
Dirac mass coefiicients B, C, E and F, derived from the complex parameters ji^, ^2, 
Ps, P2- No dependence arises from the heavy neutrino Majorana phases ^3,2 since 
they cancel out in |^?3|^ |^?2|^- 

In the limit ^2 0, where A — C — Q, but with B of order unity, the dominant 
contribution to comes from higher order terms. We can obtain it from 

(/i^g)'' - F* |^3|' (1 + \F\^) + F* (1 + F*E) Q2Ql + F* (1 + F*Ey g*gs , 
iKsf' = S (1 + F*E) Q2QI + S (1 + \B\' + \E\^) \q2\' . (4.24) 



Note that the leading term, proportional to \q^\ \q2\ ,i?> real, and in fact we did 
not ha 
terms, 



2 

not have any \B\ contributions at that order above. Hence, we consider the next 



iK^f" {K^f iKsf' « (1 + F*E) F* {k,E + K2F) Q2Ql 

+ (1 + F*Ey F {kiF* + K2E*) Q*Qs , (4.25) 
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where we defined the real parameters 

Ki = {l + \Bf + \E\^) (1 + , 



K2 = \l + E*Ff . (4.26) 



Note again that the two terms in Eq. (|4.25|) are exactly conjugate to each other 



ior E = F when the two heavy eigenstates are nearly aligned. In this limit tan ^12 
becomes maximal. Therefore, if B gives the dominant contribution, the Dirac type 
CP violation is suppressed for maximal solar angle. The CP invariant vanishes as 
well if -B = as the system effectively reduces to two generations and sin 6*13 = 
(recall that we are already in the limit A = C = 0). We then obtain 

lm[{h:^Y'{h:^f'{h:^f'] = \Bf\g,f\gsf{^^-^2) lm{Q) (4.27) 

with 

= (1 + EE*) E*(E-E) — , 

which yields 

Je \Bf{Ki~K2) Im(fi). (4.28) 



Comparison with Eq. ( [4.20]) shows then that in this case the standard Dirac phase 6 is 



a complicated function of the phases of /is, ps, p2 in the leptonic Dirac mass matrices, 
the difference between two of the Majorana phases A032 and neutrino masses. It is 
suppressed by the ratio |f?2| / l^isl, as is sin 6^13. 

Whenever only few of the parameters in the Dirac neutrino mass matrix mat- 
ter, we expect correlations between the lightest eigenvalue, the mixing angles and 
the maximal value for J^. In Appendix [A.2| , we consider the simple case where B 



dominates and the lightest eigenvalue mi vanishes; then all the observables are only 
function of B, E, F, Q2/ Q'i and we show relations among them. In this specific case, 
even allowing for the uncertainty on the phases, upper bounds can be obtained for 
sin6'i3,mee and J^. In the more general case, subleading terms and other parameters 
become important and relax any such bounds. 

4.4 Leptogenesis 

The out-of-equilibrium decays of heavy Majorana neutrinos is a natural source of 
the cosmological matter-antimatter asymmetry [33]. In recent years this leptogenesis 
mechanism has been studied in great detail. The main ingredients are CP asymmetry 
and washout processes, which depend on neutrino masses and mixings. 

It is convenient to work with a diagonal and real matrix for the right-handed 
neutrinos, which is obtained from by the phase transformation 



M 



diag (e"^"^^ , e''^^ , e-'^"") . (4.29) 
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For hierarchical heavy neutrinos the generated baryon asymmetry is dominated by 
decays of the lightest state Ni. In supersymmetric models the corresponding CP 
asymmetry is [34] 

3 Im(M?,.)Mi ^ ^n-m^* 

where the matrix elements are given, analogously to Eq. (|3.10|) , 

M., = e^^^^^p, (1 - b^) , h ^ P^^l+J'^^ . (4.31) 

PiM 

The terms involving one index 1 simplify for pi = as 

Mn = P?|l-S 

M\ 

\ M pjM ' 



The result then reads 

^1 

c 

where 

rjj = — Im 



3 Ml ( Ml\ '^ ^ p] 



1 _ ^ y J^r]. , (4.33) 



i=2,3 



M4 PjM4 + p*Mj 



(4.34) 



M PjM 

Since p|M3/(p|M2) ~ 0.2, the CP asymmetry is dominated by the intermedi- 
ate state N^, i.e., ei ~ 3/(87r)MiA/A??7|^/t>^. In any case, the phases involved, 
A013, A012 and the phases of P3,P2, are completely independent of the low-energy 
CP violating phase in the quark sector and also not so directly connected to that 
in neutrino oscillations (even if they can contribute to it). For Mi ~ 10^'' GeV, one 
obtains ei ~ 10~^, with a baryogenesis temperature Tg ~ Mi ~ 10^° GeV. These 
are typical parameters of thermal leptogenesis [35,36]. 

The strength of the washout processes crucially depends on the effective neutrino 
mass 

It 

With the efficiency factor [37] 

s/O.OleVy-^ ^^_2 



~ 10"^ ( ) ~ 10"^ , (4.36) 
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one obtains for the baryon asymmetry 

r]B ~ IQ-^EKf ~ lO-^Kf ~ 10"^° , (4.37) 

consistent with observation. So for successful leptogenesis we need a non vanishing 
Pi, Qi and in particular gi ~ §2- In such case a zero neutrino eigenvalue is only 
possible due to alignment. 

In the above estimate of the baryon asymmetry we have summed over the lepton 
flavours in the final state. In general, the CP asymmetries as well as the washout pro- 
cesses depend on the lepton fiavour, which can lead to a considerable enhancement 
of the generated baryon asymmetry [38,39]. The neutrino masses Mi ~ 10^° GeV, 
fhi ~ 0.01 eV lie in the 'fully fiavoured regime' where these effects can indeed be im- 
portant [40]. Hence, depending on the CP violating phases the generated asymmetry 
may be significantly larger than the estimate ( [1.37| ). 

5. Conclusions 

We have studied in detail a specific pattern of quark and lepton mass matrices ob- 
tained from a six-dimensional GUT model compactified on an orbifold. Up quarks 
and right-handed neutrinos have diagonal 3x3 matrices with the same hierarchy 
whereas down quarks, charged leptons and Dirac neutrino mass terms are described 
by 4 X 4 matrices which have one large eigenvalue O(Mgut)- The origin of this pat- 
tern are diagonal mass terms for three ordinary quark-lepton families together with 
large mixings (9(Mgut) with a pair of SU(5) (5 + 5) plets. This vectorial fourth 
generation though is made of different split multiplets allowing for a relaxation of 
GUT relations. The six mass parameters of the model in the quark sector can be 
fixed by the up and down quark masses. This pattern of mass matrices has several 
remarkable features: The CKM matrix is correctly predicted and the electron mass 
is naturally different from the down quark mass. 

The mismatch between down and up quark mass hierarchies leads, via the seesaw 
mechanism, to three light neutrino masses with a much milder hierarchy. Left-handed 
leptons and right-handed quarks have large mixings. This leads to large neutrino 
mixings and to small CKM mixings of the left-handed down quarks in agreement 
with observation. 

Factors 0{1) of the mass matrices are unknown, and the predictive power of 
the model is therefore limited. The neutrino mixings sin ^23 ~ 1 and sin ^13 < 0.1 
are naturally accommodated. The corresponding neutrino masses are mi < m2 ~ 
V^'^soi < "^3 ~ A/Amf^ and \mee\ ~ V^^so\ ^ O-O^ ^V. 

The elements of the mass matrices arise from a large number of different oper- 
ators. Hence, most of the CP violating high-energy phases are unrelated. We find 
that the measured CP violation in the quark sector can be obtained, even if the CP 
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invariant is suppressed by the alignment between the two Hghtest mass eigenstates. 
Due to the uncertainties of 0{1) factors no useful upper bound on the CP violation 
in neutrino oscillations is obtained in general. Some constraints can be given in the 
limited case where the number of dominant parameters is reduced, as it happens if 
the parameters A, C in the neutrino Dirac mass matrix are suppressed by the small- 
ness of the electron mass. It is indeed intriguing that in our setting the smallness 
of the reactor angle can be connected to the lightness of the electron. The model 
is consistent with thermal leptogenesis, with a possible enhancement of the baryon 
asymmetry by flavour effects. 

We conclude that mixings (9(Mgut) of three sequential quark-lepton families 
with vectorial split multiplets. a pair of lepton doublets and right-handed down 
quarks, can account simultaneously for small quark mixings and large neutrino mix- 
ings in the charged weak current and, correspondingly, for hierarchical quark masses 
together with almost degenerate neutrino masses. The CP phases in the quark sec- 
tor, neutrino oscillations and leptogenesis are unrelated. Quantitative predictions for 
the lightest neutrino mass mi and sin ^13 require currently unknown 0{1) factors in 
more specific GUT models. 
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A. Mass matrices 



We will discuss here the mass eigenvalues and the mixing matrices for the low energy 
theory in relation to the high energy parameters. 

Given a general matrix of the form as in Eq. ( |2.4| ), 



/ fii Jli\ 

/i2 /i2 

TTl = 

f:^ 

\Mi M2 M3 Ma J 

where fii,Jli = 0{vi^2) and Mj = O(Mgut), the matrices f/4 and V4 that single out 
the heavy state can be given as [14] 



1 









1 





hp 
Mi 



M2 



wMi+/£iMi\ 

^l'iM3+^^l3M4, 

M2 






1 



(A.l) 






Ml 

\ ^m2+m| 







M3 



Mia/m|+M| Mi>^ 
M^/^+^ M 



M2 



(A.2) 



with M = yJ2i ^i- III general V4 contains large mixings, while f/4 is approximately 
the unity matrix, up to terms O (v/M). Next, U3 and V3 = V^V^ diagonalise 



m' = U\mVi 



so both f/3 and V3 have a non-trivial 3x3 part only. In the following we will use 
the symbols f/3, V3 both for the that non-trivial upper-left corners and the full 4x4 
matrices obtained adding a row and column of zeros and a diagonal 1 to those. The 
effective mass matrix fh can be brought into the upper triangular form by a unitary 
matrix V3 ~ V3 such that 




m 



fh V3 



Jl2 aji2 
/is 



22 



With Vi = {mii,mi2,'fhis), the new basis is given by 



V3_ 

1^3 I 



62 



V2_ 
1^2! 



63-^2 

V2 



63 



ei = 62 X 63 . 



(A.3) 



Note that V3 corresponds to a large angle rotation for the right-handed quark fields. 

While JI3 and JI2 are real by construction, we have the freedom to choose any 
entry of the first row to be real. For concrete calculations, it is convenient to have 
7/ii real or even use the parameters as given in the basis (|A.3|) ; however, 7/ii vanishes 
in the limit fi2 — ^ 0, so for a general discussion, it is more appropriate to have Jli 
real. Here, we list the entries of fn with 7/ii real in a general form, 



1^3 



a/i2 



1^2 



pfil 



/il 



^3 



f^3f + \f^3 



1^ 

M 



IJL2HI l^2M2 + /i2M4 /i*M3 + /i3M4 



1^3 



M 



JI3M 



|/i2| + |/i2| 



1 
M2 



I - |2 

\ajj.2\ 



fiifis fiiMi + /iiM4 + li^Mi 



1^3 



/i2 

/^l I — 

/i2 



*fi3 \ _ /ilMi + /X1M4 
/^3/ M 



M4 

M 



/i2 



.*/f3 
/i3 



V'2 M 1^3 M 



7/^1 = 17/^1 



/iiMi + /I1M4 



i/^ir (i + i/^r)- 



(A.4) 



In particular, we find as well the simple expressions 

a/22 /U3^4 + /i3^3 



1^2 



7/^1 



/23M 

_ /i2 /i3 Ml 
/^2/i3 M 



/iiMi 
Ail M 



/i2M2 
/^2 M 



/il 



/i2 /i3 _^ /i2 /i3M3 - /i3M4 

/^2 /23 M ^2 JiJa 



(A.5) 
(A.6) 



These expressions vanish trivially in the limit /ii, /i2 and then we obtain the 
limiting case discussed in Section ||. As already discussed in Section H, /3/ii is inde- 
pendent of /i2. 

A.l Down Quarks and Charged Leptons 

Mass Eigenvalues and Eigenvectors. Now take the matrix fn as a starting point 
and compute the eigenvalues, eigenvectors and mixing matrices. For making things 
simpler, consider for the moment all the parameters as complex, even if actually /Is, 
/i2, 7/ii, or /23, /i2, /il can be chosen real absorbing the phases into V3. To compute 
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the eigenvalues, it is better to consider the hermitian matrices m^m or mmK The 
first option simply gives 



lAiip7 |/i2p + |/2iP |/l2p« + |/2ip/3 

J/2i|V* |/i2|V + |/2iP/3* |/i3p + |/i2|>P+|Pin/3'^ 



(A.7) 



Then the determinant is simply 

det (m%) = |det (m)|' = h\^\l2i\^\l22\''\jl3\' 



(A.8) 



and is only non- vanishing if ^Jli ^ 0. 

The eigenvalue equation is a cubic equation; to obtain the dominant terms, we 
expand around 7 = 0. In this case the equation reduces to a quadratic one with the 
solutions 



A2/3 = 2 0/^3r + 1^2^(1 + i«r) + + m')] 



(A.9) 



± 



So in this limit, we have eigenvalues at lowest order 

A3 = i/23r+i/22i>r+i/2iri/?r+0 



A2 = |AZ2r + |A^ir-0(^), 



Ai = . 



(A.IO) 



We can also compute the first correction to the zero eigenvalue simply as 



Ai 



det (m'^'m) 



V1IV2IV3P |7lViPl/"2p lA«i|«lA«2| 



A2A3 \M\\M' + \M') \M' + \M' 

This means that for vanishing //i we have 



mM' ■ (A.ii) 



(A.12) 



Using the eigenvalues, we can also solve for the mixing matrices at lowest order. 



/l 








1 




|p2pa* + |iCti|2^* 


l/^3| 





|p2pa+|piP/? 

IA*3| 

1 



(A.13) 



where we must recall that we had already acted on the mass matrix with a large 
angle rotation V3, so the above is just a small correction to it. 
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For the left-handed quark fields, we have instead at leading order 



£1 



/ii 

JI3 



(A.14) 



Since the up quark mass matrix is already diagonal, this last mixing matrix cor- 
responds to the CKM matrix. From UlfnV^ = m*^s, we get Vckm = U3, so for 
a = P we have the prediction Vtd = (a* — /?*) J^l/jis = at leading order, and the 
CP violation vanishes! On the other hand, Vub has the right order of magnitude as 
we thought. 

Quark Masses and Mixing Angles. We can reproduce the observed quark mass 
eigenvalues and mixing, that satisfy the relations 



m„ : nic : ^ A : A : 1 
rrid '■ nis : ~ A^ : A^ : 1 



(A.15) 



where A ~ Vus ~ 0.22 is the Cabibbo angle. In fact, if we assume 

/ii : /i2 : /is ~ A^ : A^ : 1 , 
yUi : /I2 : At3 - : A^ : 1 , 



(A.16) 



it gives correctly 

\Vus\ ~ 

moreover, 





AH 






V'l\ 


1/^1 




1/^3 



cb\ 



(A.17) 



I7l 



a\ 



/i2 


l/^ll 


iU2 


/is 



mb ~ AA nib — A nib ■ 



(A.18) 



Again Vtd is suppressed by the difference of a* — P* ~ 1^2! 1^2^ A^i/a^i? ^is is the Jarlskog 



invariant, Jq. 



Low-energy CP violation As discussed in the following Appendix, we can ex- 
press the low-energy CP violation in the quark section via an effective Jarlskog 
invariant. We calculate this invariant, using Eqs. (|A.4| ). The dominant terms are 
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displayed in Eq. 



the complete expression reads 



Ja 



JliJl'^fli M2M3 




M2 ■ " ^ 

M1M3 
M1M2 



fisUl , MsMi |/i3 



M2 



Im 



(A.19) 



+ 



+ 



M2 



1 - 



- Im 






/is 


Mi 






^ Im 






Mi 








M2 





/i3/ij M3M4|/i3| 



/i3/il 

^1 l/isi 



M2 



/is 



Im^ 

Hi 



M3 
M2 Jll 



/ii 



Ma 



/^l 



M2 i^i 

Ms 



Im 



/i2/il 



/i2/il 

/is/it 



Im :c — Im 

M /il M /is/ii 



M3M4^ fiS M2M. fl2 , M2M3^ /i3/i2 

— Im — — Im — H — — Im ^ ^ 

M2 /is M2 /i2 M2 /iS/i2 



Charged Leptons. The charged leptons show a different hierarchy than the down 
quarks, we have in fact 



rUe '■ rUfj, : rrir 
md-.nis: nib 



^ : : 1 
~ : : 1 . 



(A.20) 



The discrepancy can be solved with a smaller value for (/i2/ii)e, compared to (/i2/ii) 
As an example, we choose fi^ ~ A"^ and Jll ~ A'^~^ such that 









+ 


a 2 



I /ill 



l/iil 




1/^2 1 


i^S 



,2xS-4, 



A 



5-6 



(A.21) 



Regarding the rotations, the large V4 rotation acts now on the left-handed fields, 
but it has to act on both the charged leptons and the neutrinos, so it has not a large 
effect in the charged current. There is, however, an effect coming from the mismatch 
between the two Vs's in the charged leptons and neutrino cases. 

A. 2 Neutrinos 

The charged lepton mass matrix is eventually diagonalised via V3 = ^3^3' and f/3 as 
the down quark matrix. For the light neutrino Majorana mass matrix, given by 



v ( D\~^ ( N\-^ D 

m^ff = — ym j ym j m , 



(A.22) 



we can neglect the rotation f/3 of the right-handed fields as this transformation 
cancels out. U4 does in principle rotate the RH states, but its effect is suppressed as 
long as Mi < M. Regarding V3, we do not expect it to be the same for both charged 
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and neutral leptons, so the mismatch between the two provides flavour mixing in the 
neutrino sector. 

The neutrino Dirac mass matrix can be written after the large rotation V3 that 
bring the charged lepton mass matrix into triangular form as 



= m^ ¥3= I Bp2 Ep2 P2 



'Api Dpi pi' 
Bp2 Ep2 p2 

^Cp3 Fp3 



where 



Pi 



1 1 



Pi 



P3 M2 

1 1 

1 1 
P3 = — 



P3 M2 



P3 



P2M2 



P3 



and, using the notation Ma(, 



A 

B 

C 
D 

E 



1111 

Pi P2 P3 M 

P2P2 - P2P2 P3MI 

P2P2 P3 M 

P3P3 - P3P3 P2 Ml 
P2M 



P1H2P.3M1 - pi fl2P3M2 + pi2 ( P3M3 - /i3M4 



P3P3 
1111 

Pi P2 pi M2 



1111 
P2 P2 pI M2 



Pi 



P2 IPsI ^12 + PIP3M2 P3M3 - p*M4 



+ PlMi 



/ia/ig (P3M4 - /i3M3 ) + /i*M2 (/ig + |/i3|^) 



P2 



+ 



Jl2 l/isl ^12 + PIP3M2 (p3Ms - filM^j 
P2 P2PIM2 (jl3M3 - HsMi 



+ 11*2 [pilMl^ - 2 |/i3|/X3M3M4COS^3 + \P3\" M- 



(A.23) 



14 



1111 
P3 P2 pI M2 



(P3P3 - P3P3) [P2P3M/2 + p^Ms [pi^M^ - pilMij\ . (A.24) 



Note that we are here projecting the neutrino flavour states into the basis defined 
by the charged leptons as in Eq. (|A.3|) . So we can immediately see that if the neutrino 
flavour vectors are aligned with the charged leptons B,C,F should vanish and the 
neutrino mass matrix would become triangular as well. This corresponds to having 
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exactly the same hierarchy in the rows of the charged and neutral lepton Dirac mass 
matrices, i.e. ^ = ^. We do not expect such alignment since the parameters 

/^i pi 

Pi,Jii are generated by different operators and not related by any GUT relation, as 
can be seen from Eq. ( p.3c| ). We will consider in the following the case where the 
neutrino hierarchies are similar to those of the down quark matrix, while the charged 
leptons differ due to the lighter electron mass. Of course even more involved scenarios 
are possible. In the following we neglect as well corrections coming from the final 
diagonalisation, since the entries of are suppressed by {^2/^2)'^ ^ 0.01. 

Mass eigenvalues and eigenvectors. We need to compute the eigenvalues of 
the neutrino mass matrix and the first step is again to compute the determinant of 
the matrix m'^Q. Note that this is a symmetric matrix, but not real. Therefore the 
eigenvalues are in general complex and the matrix is diagonalised using a unitary 
matrix as 



Vj^es^y = diag (mi, m2, m^) . (A. 25) 

Consider for the moment just the absolute value of the eigenvalues and then see that 
we have the relation^ 

(det (m^)? 

detK,)^-i^A_^. (A.26) 

The last determinant is simply the product of the heavy neutrino masses, while the 
first one is given by 

det(m^) = P1P2P3 [{F -E){A-B) + {D- E){B - C)] . (A.27) 

In order to have three non-vanishing eigenvalues, we need all Pj 7^ and at least one 
of A, B, and C different from zero. Also the three vectors corresponding to the rows 
of the Dirac matrix must not be aligned with each other. So we obtain 

m,m2rm = -Q1Q2QZ [{F - E){A - B) + {D - E){B - C)f 

pi P2 P3 1 1 1 Ml 
= -QiQ2Qz — — — ^^^— (A.28) 
Pi P2 P3 P2 Pi M2 M 

X {pI IPsf + 2 IP2I P2 IpsI M2 P3M3 cos (6*2 - 63) - PsMi cos 6*2 
+ IP2I' [ps^is - 2 IpsI P3M3M4 cos ^3 + |P3|' Ml4 

for pi = 0, where Qi = c'^'-'^yf /Mi. 



^Note that for a n x n mass matrix, the minus sign on the r.h.s. gives a (—1)" contribution. 
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Singling out the heaviest mass eigenstate. In the case when ^3 ^ ^2,1, it is 
easy to single out the heaviest eigenstate: 

1 



\T 



(<3) 



yi + |F|2 + |C|2 

and the mass eigenvalue to lowest order is given by 



(C*,F*,1), 



rrin 



-g3{l + \Ff + \Cf). 



;a.29) 



(A.30) 



Then up to a rotation in the 12 submatrix, at lowest order the mixing matrix can be 
written as 



-CF* 





1 



\/l+\F\^-\C\- 



A 



-C -F 1 



(A.31) 



this is the basis which gives decoupling of the first eigenstate in the limit of vanishing 
C. From this matrix, we can directly read off the dominant part of the mixing angles 
with the heavy eigenstate, ^23 and ^13. The charged lepton mass matrix is nearly 
diagonal, so we can actually relate with good accuracy the first row to the electron 
neutrino flavour. The left-handed charged lepton flavour eigenstates are given as a 
function of the mass eigenstates by 



and therefore the neutrino flavour eigenstates correspond to 



t 



(A.32) 



(A.33) 



where V3 cancels out as it acts equally on the whole lepton doublet; moreover, as we 
have seen, is limited to the 23 corner and does not modify the electron entry. We 
use here the convention of [32], and deflne the PMNS matrix as 



C13C12 

-S12C23 - Ci2S23Sl3e"^ 
i5 

S12S23 ~ Ci2C23Si3e 

'10 ON / ci3 
C23 S23 1 I 

,0 -S23 C23/ \-s13e*' 




C13S12 Si3 
C12C23 - Sl2S23Sl3e"^ Ci3S23e 
-C12S23 - Sl2C23Sl3e*'' Ci3C23e 

si3e-^n / C12 S12 0\ /I 
10 -S12 C12 e^«2/2 

ci3 / V 1/ \0 e 



1 

e^^2/2 
e^^3/2 , 







i(5+«3/2) 



(A.34) 



where 
phases. 



= cos9ij, Sij — sin9ij, 5 is the Dirac phase and ^1,2 are the Majorana 
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So we have at lowest order for ^^13 that 

(K°)i3 = Sin ^^13 - , . (A.35) 

1 + + |C|' 

This gives us directly a constraint on the parameter C from the upper bound on 
|sin^i3| < 0.1: 



\C\ ~ + + |C|'|sin^i3| < 0.1-^1 + (A.36) 

Then since the mixing with the first flavour is small, the atmospheric mixing matrix 
is given simply by requiring the 23 corner of the matrix in Eq. ( |A.31j ) to give 

_ / cos ^23 sin^23e-*«^«\ , 



^— sin 6'23e*^23 q^^ 

So considering the 23 sector, we get, again at lowest order, 

^23 = arg (F) , 
tan ^23 = \F\ . (A.38) 

To have large mixing angle tan26'23 > 3 [29,31], we must restrict |F| between 

0.7 < |F| < 1.4 . (A.39) 

Such a value is natural in the case where p3, p3 and /i3, /I3 are of the same order 
but not exactly equal, while /i2 is small. Note that even a phase difference can be 
important. Assuming simply ^ = e*'^^^ and degenerate Mi gives 

^ 2V2(l-cos.3) _ 

3 — cos uj^ 

so we obtain |F| = 1 for the maximal phase difference = it, while |F| > 0.7 arises 
in the wide interval 0.26 vr < u;3 < 1.73 vr. Hence, a nearly maximal atmospheric 
angle is natural even for the most simple choice of parameters. Of course, more 
solutions are possible for the general case. 

Thus in order to reproduce the observed pattern of mixing parameters, C has to 
be small, while \F\ is nearly unity. We can use the maximal value for \F\ and the 
experimental bound on 613 to derive an upper limit on |C|, 

|C|<0.17, (A.41) 

in agreement e.g. with the ratio ^ necessary to have a small electron mass. Note, 
however, that we can obtain signiflcant corrections from f)2,i 7^ 0. 
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Light eigenstates and solar mixing angle. The other two eigenvalues and the 
correction to the heavy mass can be obtained from the trace and determinant of the 
matrix ("^^g)^ m^g, which can be computed in any basis. Expanding both the mass 
matrix and the eigenvalues to first order, 



we have then 
5m^ = 



tr \m\„m 



^3 — ^3 + <^"T.3 while mi,2 = 5rni^2 , 
] 



(A.42) 



S3" "81:2 



\mi\ + \m2\ + {Sm^l —tr 



gi,2^ ' 



7712 



|detKff)r 



(A.43) 



Choosing the basis appropriately, the relations can be simplified to give 



m 



1^/33 ' 



mi + 777.2 



7772 P^l 



tr 



331 



\QiQ2 



^2\{F-E){A-B) + {D-E){B-C)f 



{i + \Ff + \c\y 

We will give the result of these expressions for vanishing C and Qi — qg2: 

{1- FEf + q{l- FDf 



(A.44) 



5mz ^Q2- 



tr 



|"72| \mi\ ~|^2| \q 



+2 \BE + qADf + 2 |S + qA\^ + 2\E + qDf] , 
4, ,2\A{F-E) + B{D-F)\'' 



Then the mass splitting which should generate the solar oscillations is given by 

5ml^^ = \J {\mi\^ + |7772|^)^ - 4 |7n2|^ 



(A.45) 



(1+ |F| 



{ [(1 + (|1 + g|2 + + gL,2|2 ^ 1^2 ^ ^^2|2 

+ 2 + g^Dl' + 2 |5 + + 2 |£; + gL>|') 
- |(l-F£;f + g(l-FL')2|'l' 

9 ~l 1/2 

-4|g|'(l + |F|=^) |^(F-^) + 5(L>-F)|'| . (A.46) 
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So the solar neutrino mass splitting can be matched even in the case g = or 
A {F — E) + B {D — F) = 0, i.e., when the lightest neutrino is massless. However, 
we do not expect the first limit to be realised, if we assume the same hierarchies 
between pi as in the jji in the down quark sector, while for Mj as the up quark 
sector. In that case we have in fact \g2\ ~ |^?i| and the two lighter eigenvalues are 
similar in scale, mi ~ m2 — a/ ^"^soi- On the other hand, the determinant could be 
suppressed by alignment, i.e., for \A {F — E) + B {D — F) \ ^1, and could give us a 
hierarchy also between the two light eigenvalues. 

We can then compute the solar mixing angle and the first order corrections to 
the Ks mixing parameter. After rotating with the matrix, we can estimate the 
solar angle by using only the 12 part of the mass matrix; for C ~ the matrix is 
given by 



Taking the solar mixing matrix as in Eq. (|A.37|) with 6*23,^23 — * ^12,^12 we obtain 



^ ("^^1,2)12(^^1.2)11 + ("-^^1.2)22(^^1.2)^2 
1(^^1,2)12(^^1.2)11 + ("^^1.2) 22(^^1.2) 1:2 I 



tan 26* 



12 



2 


("^£.1,2 




n + 


("^51,2)22 


^"^^1,2)12! 




("^^1,2)22 


2 


("^^1,2)11 


2 



2^^/l + \F\'\A^\ 



where, for q = gi/g2, 



M = [B{E -F) + qA{D - F)] [B^ + qA^ (l + 

+ [{E -Ff + q{D- Ff] [B {E - F) + qA {D - F)]* , 

V = \{E - Ff + q{D - - \B^ + qA^\^ (l + . 

In order to have a large solar mixing angle, either Aq ov B must not be small 
compared to E — F and D — F. But since A, C oc we are led to the case 

A = C~0, B = ^^ = 0(1). (A.48) 

P2 M 

Then we can neglect the terms proportional to A and we have simply 

, / —-\Bfil + \Ff) + \E-Ff + q{D-FY^-^^ 

tan2^i2 = 2 LB E - F Jl + Ff^—- ^ ^-^^ . 

^ \{E-Ff + q{D~Fff-\B\'{l + \FfY 

(A.49) 
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This formula simplifies further if we neglect the q{D — F) terms as well.'' Then using 
general trigonometric formulae leads to the expression in Eq. (|4.11|) , 



tan ^ J^^p\ v/l + l^l' • (A.50) 

Taking the experimental value for the solar angle, tan^ 6*12 = 0.45 ± 0.05, gives us for 
|F| ~ 1 the range \B\ ~ (0.45 - 0.50) \E - F\. 

We can also compute the corrections of order ^12 to the other two mixing angles, 
that we have discussed in the lowest order. In fact, since /i2 <^ At2, the contribution 
( IA.35D is small and the leading contribution to 6*13 comes from the Bq2 term, 



(r«),3 = sin^,3 - ^^^^M -\B\^- 0.2 \B\ . (A.51) 



{^ + \F\7 



\Q3\ 



So even for vanishing leading order, we expect the first order term to bring ^13 near 
to the experimental bound. Note that it is the large solar angle that naturally gives 
^13 ~ Q2/ Qs'i in our model it seems pretty difficult to suppress this angle to much 
smaller values, apart if there is a tuned cancellation between zero and first order. 

The corrections to the atmospheric angle are of the same order f)2/f?3 and do 
not have a large effect since we need in any case large parameters in the 23 sector. 
This small shift can in fact be easily compensated by a small change in the value of 
F, especially since we do not have any particular symmetry in the model imposing 
F=l. 



Sum Rules for B dominance and vanishing mi . We have seen in the previous 
paragraph that in case of vanishing C, A and gi, simple expressions can be obtained 
for all observables as functions of only few parameters B, E, F and ^3,2- Then it is 
possible to obtain relations between the different observables. 



tan ^23 = \F\ 



\B\ 



tangi2= + 



sin 6 



13 



\B{EF + 1)\ |g2| 
(l + \Fff'\(^^\ 



"Note that taking A ^ C = D — F = gives a zero determinant for the neutrino mass matrix, 
so this case applies when the lightest eigenvalue is suppressed compared to the solar mass scale. 
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Now we can write the following relation, 



sin6'i3 5m, 



''atm 



\E-F\ \EF + l\ 



tan ^12 5msoi 



[(1 + (1 + + |E-F|^tan2^i2] - \l - EF\ 



(A.53) 

To estimate its value, we can use the fact that |F| ~ 1 and vary only \E\ and the 
phases of -E, F. We obtain then a maximal value of the r.h.s. for EE = 1 so that 



sin 6^13 < 



5m. 



sol 



tan^^ 



12 



0.09 . 



(A.54) 



Of course, the angle 6*13 can always be reduced by an appropriate choice of the phases 
and in particular for E = F, so that there is no lower bound in this type of models. 

The effective neutrino Majorana matrix, which is relevant for neutrinoless double 
beta decay, simplifies such that 



1^11^2 

5mso\ 



tan^ 9,2 \E-F\ 



+ (l + l^l^) +tan'^9i2\E - FfY - \1 - EE 



(A.55) 



Again varying the phases and the modulus of E, we find the maximal value for 
EE = -1, 



\mpp\ < 5m. 



tan^^ 



12 



sol 



0.43 5m. 



sol • 



y'2 + tan^ 612 

Moreover, we can give a simple relation between mee and the reactor angle, 

\mpA \E - E\ 



(A.56) 



5m 



atm 



+ 11 



sin 6*13 tan 6* 



12 



(A.57) 



Note that the singular value for EE + 1 = corresponds to a vanishing reactor angle. 

We can even derive a maximal value for the Dirac CP violation for this case. 
From Eqs. (|4.21j ) and (|4.27|) we get 





B 


^ (Ki — K,2) (Q) 


(1 + |F 
\E-E 


4 


a + 


E 


')2 (^i + \E\^ + \B\^Y EEf 
tan^ 6*12 (1 + tan^ 6112) Im {Q) 



(A.58) 



1 + \Fr [{1 + \Ef) (1 + \Ef) + \E- Ef tan2 6,2]- |1 - EEf 
5m,oi \E - tan2 6^2 (1 + tan^ 6^2) Im [(1 + EE*) E* {E - E) e^^^a] 



5m 



atm 



1+ li^l 



((1 + |F|^) (1 + |E|^) + |E-F|\an2^i2) -\l-EE\ 



3/2' 
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where A23 is the phase of Q2/ Q?,- Again, the prefactor is maximal for EF — —1 and 
E = —F, giving 

Ji <^^ ^ sinA23 < 0.06 . A.59 

Sm^tm 2 tan ^12 (2 + tan^ OuF 

Here the imaginary part is only given by the phase A23, but in more general cases the 
phases of E and F will play a role as well. So even for the CP violation in the leptonic 
sector, the model displays a suppression given by the ratio of the mass eigenvalues. 
Contrary to the quark case, however, the CP violation is not proportional to the 
smallest mass eigenvalue, but it can be non-vanishing even for mi = 0. 

B. CP Violation and Weak Basis Invariants 

For completeness we discuss here the CP invariants in the case of an additional 
vectorial state. We prove that if the additional state is much heavier than the 
electroweak scale, the low energy CP violation can be expressed by the Jarlskog 
invariant defined from an effective 3x3 down quark mass matrix. 

The transformation of a Dirac spinor '^(t, x) under parity and charge conjugation 
is given by 

F^it,x) p-i = 77p7V(t,-f), 

Cij{t,x) c-'^'ncCip{t,xy, 

where rjp^c are non-observable phases. The matrix C obeys the relation 7^,^(7 = 
— C7J. Since the Lagrangian is a Lorentz scalar, it only depends on fermionic field 
bilinears. Thus, we deduce the CP transformation for such terms, 

CP v^iV'i (CP)"' = , 

CP (CP)-' ^ -^j^'A : 

CP AYA (CP)"' = -Ali^A , 

CP AY7% (CP)"' = -A7>.7'A ■ 

Note that the operator d/^ transforms under CP as ^ d^. 

Quark Sector. In the Standard Model, it is easy to verify the existence of the CP 
symmetry in the Lagrangian, up to mass terms. In general, the quark mass terms are 
CP invariant if and only if it is possible to find a weak basis transformation which 
realises 

= WlH^'wI , H'^* = WlH'^wI , (B.3) 

where = M"'"^ (M"''=')^ It follows that 

Wl [H^, Ha] Wl^- [H^, HaV , (B.4) 
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such that, for r odd, 



tT[H^,HdY = (B.5) 

is a necessary and sufficient condition for CP invariance [42]. 

The case of r = 1 is trivial: the trace of a commutator [Hu, Ha] is zero. For 
r = 3 and three generations, we have 

JsM - tr [H^, H,f = 6. K - -c) (-? - <) i^l - ml) 

/ 2 2\ / 2 2\ / 2 2\ T K^-^J 

where the quantity Jg does not depend of the mass spectrum, and can be related, up 
to a sign, with the CKM matrix, V, as \ Jq\ = |Im(Vi2l^i*3V"2*2V23)|. We conclude that in 
order to have CP violation, we need to have Jg ^ 0. This quantity is the lowest weak 
basis invariant which measure CP violating effects and it has mass-dimension twelve. 
Apart from CP violation in the strong interactions, there is no other mechanism in 
the SM which can generate CP violating effects if Jg = 0. Note that in the chiral 
limit, rriu = rrid = = 0, we do not generate CP violation even if Jg 7^ 0. 

In the literature, the lowest weak basis invariant is called Jarlskog determinant 

[28], 

det [H^, H,] = 2^ {m\ - m^) (m? - m^) {m\ - ml) 

12 2\ I 2 2\ I 2 2\ T V 

which is equivalent to the Eq. ( |B.6| ).'^ The Jarlskog determinant is only applicable to 
the case of three generations, in contrast to the more general invariant in Eq. ([B.5|). 



Now let us add a down quark isosinglet. The gauge couplings to quarks and 
their mass terms are = 1,2,3 and a = 1,2,3,4): 

= --^ {uurdLi W;l + h.c.) - eJ^^A^ 

- - — ^^-7— {uLil^ULi - din'^du - 2 sin^ 6'w ^em) (B.8a) 
z cos fw 

^M = - {uu M'J URj + dLi M'f dn^ + dL4 ml dna) + h.c. (B.8b) 

where the matrices M^, and are of dimension 3x3, 3x4 and 1x4, respectively. 
The electromagnetic current is given by J^j^,j = ^u'y^u — ^d'j'^d. 

The most general weak basis transformation consistent with the Lagrangian of 
Eq. (|R8D is: 



^For any 3x3 traceless Hermitian matrix M one has: tr Af'^ = 3 |Af | 
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where Ul and are 3x3 unitary matrices, while f/^ is 4 x 4. Once we diagonahse 
the mass terms, the Lagrangian reads 



•^w = [uLiY (V^ckm)^^ dia W+ + h.c] - eJ^^A^ 



9 



ULiYuLi - diaY [ VcKM^CKM ) C?L/3 " 2 sin^ 6'w 

af3 



M 



2 COS Oy^ 

[uLi Dui URi + die Dda dna) + h.c. 



EM 



(B.IO) 



where Vckm = U^^Uf is a 3 x 4 matrix. The number of independent phases which 
are related to CP violation is, for = 3, 

ncp = (iV + 1) - ^A^ (A^ - 1) - 2iV = iA^(A^ - 1) = 3 . (B.ll) 

With the matrices as defined in Eq. (p^SbD and = M^Ml, = MaM^, and 
hci = Mdm\, we can write down a set of weak basis invariants. 



/i = Im tr HuHdhdh\ , 

h = \mtiHl[H^,Hd] hdh\, 

h = \miiHlHlhdh\, 



h = lmiiHlHdhdh\ , 
U = ImtT HuH'^hdhl , 
h = lmtTHl[H^,Hj] hdh\ 



d ' 



h = \mtiHlHdH^H^ 



d 1 



(B.12) 



representing a set of necessary and sufficient conditions for having CP invariance in 
the quark sector [43]. 

In our model, Hd and hd read 



Hd 



I |2 I ~2 



( fiiMi + JIiMA 
H2M2 + /I2M4 



(B.13) 



Since Hu and Hd are real, Ii vanishes. The remaining invariants are in general 
different from zero; the dominant terms are 



h = -ml {^l + /i^) /i3M4 Im/ig 



I2 = mf h 



-m. 



h = -ml {fil + ul) (/ii + /i2 + + f4) l^aMsM^ Im /xg , 
/g = -ml {Jll + Jll) {Jll + Jil + Jil + l4) Ji^M^Mi Im /xg . 



(B.14) 



Hence, CP is generally violated even by the presence of a single complex parameter 
/is. Note that this case is not equivalent to the chiral limit because both the charm 
and strange masses are different from zero, rric cx ^2 and ~ /I2 (albeit ^2 
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JI2). As we might expect, the invariants vanish if all quarks of the first and second 
generation are massless. 

Now we single out the heavy eigenstate with the rotations V4, U4. While the 
action of V4 leaves the invariants unaffected, strongly modifies them and reshuffles 
terms from one to the other. In fact after this transformation, hd vanishes to lowest 
order and survives only at order O^vj^^r/M^); then in the new basis all the invariants 
involving hd, i.e., Ji — Iq are suppressed by u|^/M^ and vanish for M 00. On the 
other hand Ij is now non- vanishing and given by 

I', = lv^iTHlHfH^{Hff , (B.15) 

where H'f = mm) (see Eq. ( p.7|) ). Note that also changes the weak interactions, 

5^w = ua^ {Ui - 1),4 iy+ + diT'^ (ulUi - 1) + h.c, (B.16) 

V 2 v / j4 

so we expect both CP violation and CKM unitarity violation from these terms as 
well. However, the mass of the heavy state is O (Mqut) so that the contributions to 
low-energy processes are suppressed by a factor Mew/^gut and are negligible. 
Hence, at the electroweak scale, we are left to consider the single invariant 

I', = \mtiHlHfH^Hf\ (B.17) 

which corresponds to the usual Jarlskog invariant Jq for three generations, but com- 
puted for the effective quark mass m. 

Lepton Sector. As discussed above, we can ignore the heavy states for low-energy 
CP violation and use the effective 3x3 Yukawa matrices instead. 

In the SM, extended by right-handed neutrinos, we have three mass terms for 
the leptons, 

■^M = - (^^Li m'J CRj + ULi "^D ^Rj + ^ ^RiC m'^ i^Rj^ + h.c. (B.18) 

In analogy to the quark sector, invariance of the mass terms under CP transformation 
requires 

U^meV = ml , U^moW = , W^^m^W = -M^ , (B.19) 

where U, V, and W are unitary matrices acting in flavour space. Deflning h = m^^mn 
and H = rn^j^m^, we obtain 

W^hW = h\ W^HW = H*. (B.20) 

Now we can write down the weak basis invariants 

jf = Imti h H m*j^h*mN, I2 = lva.tr h H'^m*j^h*m]\[, 

4 = Im tr /i H^m%h*mNH; (B.21) 
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for the three further invariants, substitute h = rri^^me'm\mi:) for h [26]. In the basis 
where the right-handed neutrino mass is diagonal, one obtains 



/{ = M1M2 {Ml - Ml) Im h 

+ M1M3 [Ml - Ml) Im hl^ + M2M3 [Ml - Ml) Im h. 
li = M1M2 {M^ - Mfl Im h 

+ M1M3 (M^ - ) Im hj^ + M2M3 {Ml - M^) Im h 
ll = Ml Ml {Ml - Ml) Im hj^ 

+ MlMl {Ml - Ml) Im hj.^ + MImI {mI - Ml) Im hj^ 



2 

12 

2 

23 ' 

2 
12 

2 

23 5 



If none of the Mi vanish and there is no degeneracy, the vanishing of Ji, I2, and Is 
implies the vanishing of Im/ifgi Ii^^i3i and Imh^^ for CP invariance. 

Note that in our model, ttid stands for the effective 3x3 part of the Dirac 
neutrino mass matrix, m/), as given in Eq. ( [A.23| ). Then we obtain from Eq. 



h2 = A* Dpi + B*Epl + C*Fpl , 
hi3 = Apl + Bpl + Cpl , 

/i23 = D'^pI + E*-pl + F*pI . (B.22) 

The coefficients A, . . . , F are displayed in Eqs. (|A.24| ). They are generically complex, 
so we do not expect CP to be conserved. 

As in the quark sector, these invariants are rather general and give the necessary 
conditions for the presence of CP violation. On the other hand, only few of the phases 
remain important also in the low-energy limit. In our case, to study the low-energy 
Dirac invariant, we can use the analogue of the Jarlskog invariant, 

as discussed in Section ^73| . Here, /igg- = {m^^)^ m^^ and and are the 

products of the mass squared differences of the light neutrinos and charged leptons, 
respectively. 
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